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We present theoretical and numerical evidence to show that self-defocusing nonlinear optical propagation can be used to 
compute Euler fluid dynamics and possibly Navier-Stokes fluid dynamics. In particular, the formation of twin vortices and the 
Karman vortex street behind an obstacle, two well-known viscous fluid phenomena, is numerically demonstrated using the 
nonlinear Schrodinger equation. © 2008 Optical Society of America 

OCIS codes: (200.0200) Optical computing; (190.5530) Pulse propagation and solitons 



1. Introduction 

A. Philosophy of metaphoric computing 

Nonlinear dynamical systems, such as weather, plasma, and 
the economy, are ubiquitous in nature and everyday's life, yet 
such systems are typified by their highly complex and chaotic 
behaviors, making them notoriously difficult to study theo- 
retically, experimentally, and numerically. Analytic solutions 
of nonlinear systems are rare, experiments are often too in- 
flexible or impractical, and numerical simulations must take 
into account a large number of data points in multiple dimen- 
sions in order to accurately model a problem of interest, such 
that even the fastest supercomputers today would take days or 
weeks to simulate relatively simple nonlinear dynamics that 
a physical system exhibits in seconds. 

On the opposite side of the same coin, we can regard the 
physical system as a computational device that computes its 
own dynamics at a speed unimaginable by supercomputers. 
The key to harnessing this tremendous computing power of a 
physical system is therefore to make it compute other inter- 
esting problems of the same order of complexity. Of course, 
a conventional digital computer is itself a physical system, 
but it makes use of complex semiconductor physics to com- 
pute elementary logic operations, and in doing so, discards 
a large amount of information that is considered extraneous. 
In this perspective, a digital computer is an extremely inef- 
ficient computing device, as it only utilizes an exceedingly 
small amount of the full computing capability potentially of- 
fered by its physics. The advantage in this case is the flexibil- 
ity in cascading different logic operations for general-purpose 
computing, but as evidenced by the difficulties in the numer- 
ical simulations of nonlinear dynamical systems, this ineffi- 
cient computing method is often inadequate. 

In order to make full use of the computing capability of- 
fered by a physical system, we hereby propose the concept 
of metaphoric computing, which makes use of a more ex- 
perimentally accessible nonlinear dynamical system to simu- 
late another nonlinear dynamical system. An example of this 
computing method is a wind tunnel, in which a small-scale 
fluid experiment is performed to simulate large-scale fluid 
dynamics, by virtue of the scaling laws inherent in fluid dy- 
namics. Metaphoric computing, however, is not restricted to 
the use of similar physical systems to simulate each other. 
In this paper, we show in particular that nonhnear optics 



can compute fluid dynamics as well. An optical beam in- 
herently holds three-dimensional spatiotemporal information, 
and nonlinear optical propagation computes the evolution of 
this large amount of information simultaneously at the speed 
of light, promising substantial parallelism and speed for com- 
puting. Although the use of nonlinear optics for digital com- 
puting has not been as successful as the use of solid-state elec- 
tronics, forcing optical beams to compute binary logic wastes 
most of the spatiotemporal information that can be manip- 
ulated in optical beams. Instead of fitting a square peg in a 
round hole, using optics to simulate other nonlinear dynami- 
cal systems provides a natural way of making full use of the 
computing capacity offered by a nonlinear optical system. 

Fluid dynamics, the foundation of a wide variety of im- 
portant research fields including meteorology, aeronautics, 
plasma physics, superfluids, and Bose-Einstein condensates, 
is an ideal problem to solve by metaphoric computing. In- 
tractable theoretical analysis and inflexible experiments com- 
pel the use of numerical simulations, the difficulty of which 
nonetheless gives rise to a whole new field, computational 
fluid dynamics, in itself. The main difficulty is due to the in- 
herent complexity of a fluid dynamics problem, which is non- 
linear and continuously generates finer structures as the fluid 
dynamics evolves. For problems that are of practical inter- 
est, such fine structures are often orders-of-magnitude smaller 
than the size of the objects under consideration, thus requir- 
ing a large number of data points in each of the three spatial 
dimensions to be manipulated at each time step, which must 
also be correspondingly small to avoid numerical instabilities. 
An alternative method of simulating complex fluid dynamics 
that combines the speed of a fluid experiment and the flexi- 
bility of a numerical analysis is hence of great practical im- 
portance. In this paper, we show that, via a suitable transfor- 
mation, nonlinear optical propagation can be utilized to simu- 
late Euler fluid dynamics, which is known to be computation- 
ally expensive and unstable to solve numerically . We also 
provide strong evidence that nonlinear optics can simulate 
high-Reynolds-numberNavier-Stokes fluid dynamics as well, 
which include a large class of important and computationally 
difficult problems, such as turbulence. With the speed, paral- 
lelism, and configurability of optics, an "optical wind tunnel" 
may one day become a viable alternative to experiments and 
numerical analysis in the study of fluid dynamics. 



B. Correspondence between nonlinear optics and fluid dy- 



The analogy between nonlinear optics and fluid dynamics has 
been noted by many authors.'"'** Wagner et al. first suggested 
that the nonlinear propagation equation of an optical beam 
can be recast into equations that resemble the continuity equa- 
tion and the Bernoulli equation in irrotational fluid dynam- 
ics.' Coullet et al. first coined the term "optical vortices," 
which shows the analogy between phase singularities in op- 
tics and fluid vortices.^ Brambilla et al. noted that laser equa- 
tions can be transformed to a hydrodynamic form.^ Arecchi 
et al. first experimentally demonstrated the dynamics of opti- 
cal vortices in nonlinear optics.'* Akhmanov et al. called the 
rich nonlinear dynamics observed in a nonlinear resonator 
"optical turbulence."^ Swartzlander and Law observed opti- 
cal vortex solitons created via the instability of dark soliton 
stripes analogous to the Kelvin-Helmholtz instability in fluid 
dynamics.'' Staliunas showed that a laser can be described 
by the Ginzburg-Landau equation, which can be transformed 
into equations resembling the Navier-Stokes equations that 
describe viscous fluid dynamics.^ Vaupel et al. observed vor- 
tex pair nucleation by the interference of two modes in a laser 
and claimed that it was an analog of a vortex street behind an 
obstacle in a fluid flow.^ Molina-Terriza et al. also observed 
optical vortex streets in walking second-harmonic genera- 
tion.^ Roux"^ and Rozas et al.^^ studied the interactions be- 
tween optical vortices and found that their interactions resem- 
ble those of fluid vortices. Rozas et al. then experimentally 
demonstrated the fluidlike motion of a pair of optical vor- 
tices.'^ Michinel et aO^ and Paz-Alonso et al.^^ found that 
optical propagation in a cubic-quintic nonlinear medium re- 
sembles a liquid drop, and optical vortices in such a medium 
also have fluidlike motions.'^ On the other hand, nonlinear 
optics has been compared with superfluids and Bose-Einistein 
condensates, as they can all be described, to varying de- 
grees, by the nonlinear Schrodinger equation, '^'^^ commonly 
known as the Gross-Pitaevskii equation in the field of super- 
fluids.^' Pomeau and Rica suggested that the phenomenon 
of transition to dissipation in a superflow^^ can be observed 
in nonlinear diffraction.'^ Bolda et al. numerically demon- 
strated the same phenomenon in a nonlinear Fabry-Perot cav- 
ity.'^ Chiao also found that photons in such a cavity should 
obey the Bogoliubov dispersion relation for a superfluid. '^ 

The abundant amount of prior work credited above pro- 
vides ample evidence that nonlinear optics resembles fluid 
dynamics to a certain degree. In order to use nonlinear op- 
tics as a useful and practical computational tool for fluid 
dynamics, however, simply drawing analogies between the 
two kinds of dynamics is not enough. One must be able to 
show an exact correspondence, or at the very least, an ap- 
proaching convergence between a problem in nonlinear op- 
tics and a problem in fluid dynamics, in order to produce 
any useful prediction of fluid dynamics via nonlinear optics. 
Moreover, as computers nowadays have enough capabilities 
to simulate two-dimensional fluids, the mere correspondence 
between optics and two-dimensional fluid dynamics consid- 
ered in most of the prior work would not motivate the use of 



metaphoric optical computing in preference to conventional 
digital computing. A three-dimensional fluid modeling, on 
the other hand, requires a processing capability orders-of- 
magnitude higher than that available in today's supercomput- 
ers, so metaphoric optical computing would need to compute 
such problems much more efficiently to compete with elec- 
tronic computers and the Moore's law. 

In the following sections, we shall attempt to estab- 
lish the correspondence between nonlinear optics and three- 
dimensional fluid dynamics. We shall show that, taking 
group-velocity dispersion into account, nonlinear optical dy- 
namics approaches three-dimensional inviscid Euler fluid dy- 
namics in the highly nonlinear self-defocusing regime, where 
the optical intensity represents the fluid density, the optical 
phase gradient represents the fluid velocity, the nonlinear re- 
fractive index perturbation represents pressure, the propaga- 
tion distance represents time, and the temporal dimension of 
the optical pulse represents the third dimension of the fluid. 
As Euler fluid equations often exhibit high numerical instabil- 
ities, this correspondence in itself should be useful in model- 
ing high-Reynolds-number fluid dynamics away from objects 
and boundaries. In the convergence of nonlinear Schrodinger 
equation towards the Euler equations, a "quantum pressure" 
term arising from the nonlinear Schrodinger equation plays 
the role of a small parameter. As this quantum pressure term 
plays analogous roles to viscosity in the Navier-Stokes equa- 
tions, we argue that nonlinear optics should be able to ap- 
proximate viscous Navier-Stokes fluid dynamics as well, in 
the regime where quantum pressure and viscosity both play 
the role of small parameters in the respective equations. That 
said, we do not pretend that we have established the equiva- 
lence between nonlinear optics and Navier-Stokes dynamics, 
as the similarity between quantum pressure and viscosity is 
still an open problem. 

On the practical side, in cases where ideal nonlinear op- 
tics setup is not available, we suggest a split-step method that 
pieces together different optical devices to approximate an 
ideal nonlinear optics experiment. This method is very similar 
to the method proposed to simulate quantum systems using a 
quantum computer.^-' 

It must be stressed that although we focus on simula- 
tions of classical physical systems, future quantum com- 
puters that simulate quantum systems"^ would run into the 
same problem of manipulating a large amount of multi- 
dimensional information. In the case of quantum sys- 
tems, multi-dimensional quantum information, such as a 
multi-particle multi-spatiotemporal-dimensional wavefunc- 
tion, needs to be processed in parallel. Quantum comput- 
ers can naturally parallelize the multi-particle aspect, but 
there is no obvious way of parallelizing the manipulation 
of multi-spatiotemporal-dimensional information via simple 
binary quantum logic. Perhaps a quantum metaphoric com- 
puting would then be necessary, where a more accessible 
multi-dimensional quantum system is used to simulate an- 
other quantum system. 



2. Correspondence between nonlinear optics and Euler 
fluid dynamics 

A. Madelung transformation 

We now proceed to show mathematically how the self- 
defocusing nonlinear optical propagation equation, including 
the effect of group-velocity dispersion, can be transformed to 
three-dimensional hydrodynamic equations. First, we show 
how the optics equations, in the absence of optical vortices, 
correspond to inviscid and irrotational fluid equations. This 
form of transformation is widely attributed to Madelung. ^^ 
We model the paraxial nonlinear propagation of an optical 
beam, described by the envelope function \i/{z,x,y, T), via the 
nonlinear Schrodinger equation,'^-^*' 
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where z is the propagation distance, ko = 2nno/?^ is the car- 
rier wave number, j32 is the group-velocity dispersion coef- 
ficient, T is the time coordinate in the moving frame of the 
pulse, and «2 is the nonlinear Kerr coefficient. To use the time 
coordinate as the third spatial dimension of a fluid, anomalous 
group-velocity dispersion, or P2 < 0, is required. Dispersion 
can then be regarded in equal footing as diffraction if a nor- 
malized time coordinate is defined as 
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The Madelung transformation is defined as follows, 
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such that the evolution equations for the intensity, /, and the 
phase gradient, k, are given by 
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One can already see that Eq. has the exact same form 
as the fluid continuity equation, while Eq. (|8j resembles the 
Bernoulli equation,' if one regards the intensity as the fluid 
density and the phase gradient as the fluid velocity. The non- 
linear refractive index term, kon2l, would resemble the fluid 
pressure if «2 < 0, so self-defocusing is required. The last 



term in Eq. (|8} is a peculiar term that arises from optical 
diffraction and dispersion, does not exist in classical fluid dy- 
namics, and is commonly called the "quantum pressure." 

In order to compare these equations with fluid equations 
more easily, we use the following normaUzed variables, 
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where W is the characteristic size, K is the characteristic 
phase gradient, /q is some characteristic optical intensity 
of the propagation, and a is the so-called "healing" length, 
which is the length scale at which the quantum pressure term 
has the same order of magnitude as the nonlinear term on the 
right hand side of Eq. (|8}, ^ is the Mach number, which 
measures the relative strength of fluid pressure compared 
with convection, and £? is another number that measures the 
relative strength of fluid convection compared with quantum 
pressure. The normalized equations become 
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Equation (I12t is exactly the same as the fluid continuity equa- 
tion, and in the limit of ^/£ jM -^ 0, which is the highly self- 
defocusing regime, Eq. ( I13t is the same as the hydrodynamic 
equation of motion that describes inviscid and irrotational flu- 
ids. Equations ( I12> and ( I13> also admit sound wave solutions, 
which describe travelling perturbations to the density and the 
velocity. As long as the sound waves are weak, the depen- 
dence of pressure on the density is not crucial, and the use of 
self-defocusing Kerr nonlinearity is adequate. This restricts 
the correspondence to slightly compressible barotropic fluids. 
In order to model slightly compressible fluids, the opti- 
cal beam needs to have a relatively constant intensity back- 
ground. This can be achieved approximately near the cen- 
ter of a very large beam, in a large multimode waveguide 
as a container in two spatial dimensions, or in a cubic- 
quintic nonlinear medium to provide a "surface tension" to 
the beam. '^"'^'25 

B. Vorticity 

In general, the fluid velocity vector should contain an irrota- 
tional component and a rotational component. 
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where (p is called the velocity potential, and the curl of u is 
defined as the fluid vorticity. 
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The dynamics of vorticity is arguably the cornerstone of hy- 
drodynamics.^' The inviscid fluid dynamics that includes the 
rotational effect is governed by the Euler equation. 
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where P is the pressure. For incompressible fluids, ^ << 1, 
and as long as P increases with p, the specific dependence 
of P on the fluid properties is not important. Equation ( I16t 
contains the convective term u • Vu, which can be written as 



u Vu = V ( -u u 



V| -u-u 1 +fflxu 



Vxu)xu (17) 

(18) 



One can then see that the optical BernoulU equation, Eq. (I13> . 
misses the rotational component of the convective term. In 
other words, the Madelung transformation is only able to de- 
scribe the irrotational part of the fluid motion, but not the 
more important rotational part. 

The inability of the Madelung transformation to describe 
vorticity is due to the failure of the transformation near opti- 
cal vortices, where Eq. ( I13> is ill-defined. To understand this 
problem, consider a rectilinear optical vortex in polar coordi- 
nates and neglect the third fluid dimension for now. 
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where m is an integer and is called the topological charge of 
an optical vortex. The phase gradient is then given by 
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The fluid vorticity is proportional to the curl of k. 
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which resembles the vorticity of an ideal point fluid vortex. ^^ 
The motion of these vortices, however, cannot be described 
by the Madelung transformed equations due to two problems: 
k diverges when r ^ 0, so the fluid velocity u at the center of 
a vortex is infinite, and f{r) must approach r'" in the limit of 
r ^ to maintain the continuity of y/, so the quantum pres- 
sure term, with ^/p in the denominator, is also infinite near 
the vortex center. 

To overcome these difficulties, it is necessary to consider 
the motion of the optical vortices separate from the irrota- 
tional optical flow. 

C. Optical vortex solitons and point vortices 

In a relatively constant intensity background, optical vortices 
exist as optical vortex solitons.^' ^^'^^ The optical envelope 
function i/a of a vortex soliton is given by Eq. ( I19> , where 
f{r) -^ r'" for r << a, and f{r) approaches a constant for 
r >> a, where a is the healing length and also the size of 



the dark spot of a vortex soliton. In three dimensions, a vor- 
tex soliton exists as a vortex filament. We shall hereafter con- 
sider single-charged vortex solitons with »z = ± 1 only, as they 
have the lowest energy and are the most prevalent ones aris- 
ing from an experimental situation. It is also more accurate to 
approximate continuous vorticity with only discrete vortices 
with the smallest topological charge. 

Eq. J22t suggests that an optical vortex soliton resembles 
an ideal point vortex in incompressible fluids. Indeed, the mo- 
tion of optical vortices in the highly self-defocusing limit can 
be rigorously proven to behave in the same way as point fluid 
vortices.'"' "'^**"^' If one defines the position of each vortex 
filament as x^, then the fluid velocity at each point due to 
the presence of the vortex filaments in the limit of high self- 
defocusing is given by 
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where x is the normalized three-dimensional position vector, 
mj is the topological charge of vortex j, and — V^ describes 
the irrotational flow according to Eq. ( I13> . In particular, the 
motion of each filament is given by. 
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in the leading order These equations of vortex motion are 
valid as long as the separations of the vortices are much larger 
than a. For example. Fig. Q plots the intensity, phase, and 
phase gradient of two rectilinear optical vortex solitons with 
the same charge, which should rotate around each other, and 
those of two vortices with opposite charges, which should 
drift in the same direction perpendicular to their separation. 
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Fig. 1. Intensity (left column), phase (middle column), and 
phase gradient (right column) of two optical vortex solitons 
with the same charge (top row), which should rotate in the 
same sense, and those of two vortex solitons with opposite 
charges, which should drift in a direction perpendicular to 
their separation. The phase gradient near the centers of the 
vortices is not plotted due to its divergence. 

With the vortex filaments, one can define the equivalent 



voiticity in an optical beam. 
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which can be used to approximate the continuous vorticity 
of a fluid, if the number of vortex filaments is large enough. 
In this case, to include the vorticity effect in the nonlinear 
optical dynamics, one can phenomenologically patch up the 
irrotational equation of motion, Eq. ( fTst . 
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This modification of the equation of motion can be attributed 
to the phenomenon of phase slippage,^^'^^ y^gjj known in the 
field of superfluids. The use of discrete point vortex interac- 
tions to calculate Euler fluid dynamics is also a well-known 
numerical method in computational fluid dynamics. ^^ Hence, 
to simulate Euler fluid dynamics, one can approximate both 
the rotational and irrotational components of the initial fluid 
velocity profile by the optical phase and the phase singular- 
ities in an optical beam, and the nonlinear self-defocusing 
propagation of the beam would converge to incompressible 
Euler fluid dynamics in the strongly self-defocusing regime. 
One can also borrow from the well-established numerical 
techniques'^ to determine how the distribution of optical vor- 
tices sufficiently approximates the continuous vorticity in flu- 
ids. 

D. The fluid flux representation 

So far, we have shown that optical vortex solitons behave like 
point vortices in fluids when they are far away from each 
other, and this behavior can be used to approximate Euler 
fluid dynamics. However, there is no guarantee that the vor- 
tices would remain well separated in the course of the vor- 
tex dynamics. If optical vortices behaved exactly like point 
vortices, then their velocities would diverge when they are 
close to each other This velocity divergence is well known to 
cause significant numerical instability in the use of point vor- 
tices for computational fluid dynamics. '^ Another problem is 
that in three dimensions, the self-induced velocity of a curved 
point vortex filament diverges logarithmically ^ ln(l/a) in 
the limit of a ^ 0.'^ Since the optical intensity decreases to 
zero near the center of an optical vortex, the quantum pres- 
sure term, which determines the size of the vortex dark spot, 
can no longer be ignored, and the optical vortex interactions 
should differ markedly from point vortex interactions when 
their separation is on the order of a. 

To investigate the optical vortex dynamics when they are 
close to each other, the fluid velocity is no longer an appropri- 
ate quantity to study, because it diverges near a vortex center. 
The density, on the other hand, approaches zero towards the 
center. This motivates us to define an alternative finite quan- 
tity by multiplying the velocity and the density. 
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which is the fluid flux, or the momentum density. Simple cal- 
culations show that the flux is indeed finite everywhere in 



an optical beam, including the center of an optical vortex. In 
terms of the flux, the tensor dynamical equations now read''* 
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where 7, is the /th component of J, d jdxi is the /th spatial 
derivative, and repeated indices are implicitly summed in the 
manner of Einstein's summation. These equations have the 
same form as the normalized Euler equations in the tensor 
form. 
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except the quantum pressure term in Eq. ( I29> . Hence, in the 
flux representation, we have successfully avoided the problem 
of divergent quantities. Furthermore, Eq. ( I29> . in contrast to 
Eq. ( I13> . includes the correct convective term. 

The use of momentum density in the description of non- 
linear optical dynamics is more natural and appropriate than 
the use of velocity in the Madelung transformation, as the dy- 
namics ultimately evolves according to the basic law of mo- 
mentum conservation. As we shall show next, when compar- 
ing the optical flux to the fluid flux, the dynamics of optical 
vortex solitons are much more similar to that of less singu- 
lar fluid vortices than point vortices, and the correspondence 
between nonlinear optics and Euler fluid dynamics is still jus- 
tified when a is finite. 

E. Optical vortex solitons and vortex blobs 

In light of the fluid flux representation, one should therefore 
compare the flux of an optical vortex soliton to the flux of 
a fluid vortex. In an incompressible fluid, the density is con- 
stant, so the flux is proportional to the velocity, and the flux at 
the center of a point vortex has the same singular behavior as 
the velocity. Near a vortex soliton, however, the flux is finite. 
Consider the example of a single-charged vortex soliton. The 
flux near the center is given by 
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which vanishes as r ^ 0, as opposed to the divergence of 
J '^ 1 /r at the center of a point vortex. 

Instead of comparing a vortex soliton to a point vortex, one 
should hence compare the soliton to a vortex blob,^^ which 
has finite vorticity over a finite area. The vorticity of a vortex 
blob filament is mathematically described by 

(o{x,Q=2nmj j dxjy{\x-Xj\) (33) 



where 7 is a vorticity distribution function for the filament. 
The velocity near the center of a rectilinear vortex blob and 
far away from the center is 
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so in an incompressible fluid, the fluid flux of an optical vor- 
tex soliton with size a is the same as that of a vortex blob 
with size a. See Fig. |2] for a graphical illustration. The dy- 
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Fig. 2. Sketches of velocity and flux of a vortex blob and an 
optical vortex along a line across the center, to illustrate the 
similarities between the two in terms of the flux. 

namics of a vortex blob and that of a vortex soliton are also 
extremely similar. For example, the rotation frequency Q. of 
two like-charged vortex blobs approaches a constant <=<: 1/a^ 
when their separation goes to zero. Numerical simulations of 
the nonlinear Schrodinger equation also show that the rota- 
tion frequency of two like-charged vortex solitons approaches 
a constant <=<: l/a^ and does not diverge like two point vor- 
tices. ^^ On the other hand, the self-induced velocity of a 
curved vortex blob filament is given by 

dxj _ iriib Pc_ 
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where b is the unit binomial vector of the filament and p^ is 
the radius of curvative. The self-induced velocity of an optical 
vortex soliton filament is proven to be exactly the same.^^ 
Hence, optical vortex solitons act as vortex blobs, and a large 
number of solitons can simulate Euler fluid dynamics, much 
like the popular discrete vortex blob method in computational 
fluid dynamics. ^^ 

F. Numerical evidence of correspondence between nonlin- 
ear optics and Euler fluid dynamics 

The most telling evidence of the correspondence between 
nonlinear optics and Euler fluid dynamics is perhaps the 
numerical fluid dynamics simulations using the nonhnear 



Schrodinger equation by Nore et alr'^'^'' Using the nonlin- 
ear Schrodinger equation, Nore et al. numerical demonstrated 
the Euler fluid dyanmics of a jet made of an array of counter- 
rotating vortices, which exhibit sinuous and varicose instabil- 
ities.^^ In another study, Nore et al. also demonstrated three- 
dimensional shear flows and showed that numerically solving 
nonlinear Schrodinger equation is a viable alternative to Euler 
and Navier-Stokes equations for the numerical study of shear 
flows. ^^ As nonlinear optical propagation is governed by non- 
linear Schrodinger equation. The numerical experiments by 
Nore et al. show that nonlinear optics should also be able to 
compute Euler fluid dynamics. 

3. Similarities between nonlinear Schrodinger dynamics 
and Navier-Stokes fluid dynamics 

In the previous sections, we have shown the correspondence 
between self-defocusing optical propagation and inviscid Eu- 
ler fluid dynamics via a variety of methods, including the 
Madelung transformation, the incorporation of vorticity ef- 
fect due to the "phase slip" phenomenon, the fluid flux rep- 
resentation, and the comparison between optical vortex soli- 
tons and vortex blobs. Even though viscosity plays the role of 
a small parameter in most interesting fluid dynamics prob- 
lems, its effects are of paramount importance near a "no- 
slip" boundary and in the dissipation of eddies, in which 
cases the viscous Navier-Stokes equations should be used. 
In this section we shall present evidence that the nonlinear 
Schrodinger equation exhibits many of the same behaviors 
of viscous Navier-Stokes fluid dynamics, and in each case, 
quantum pressure plays an analogous role to viscosity. 

The normalized Navier-Stokes equations in the flux repre- 
sentation are given by 



(37) 
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where the last term in Eq. (ySjl is the viscosity term and Si 
is called the Reynolds number, which describes the relative 
strength of convection compared to viscosity. 






(39) 



where U is the characteristic velocity of the fluid system, L is 
the characteristic length, and v is the kinematic viscosity of 
the fluid. Comparing the viscosity term in Eq. ( I38> with the 
quantum pressure term in Eq. ( I29> via a dimensional analy- 
sis would suggest that an analogous optical Reynolds number 
would be defined as 



S^ = KW, 



(40) 



where, to recall, K is the characteristic optical phase gradi- 
ent, and W is the characteristic size of the optical experiment 
setup. The optical Reynolds number thus roughly measures 



the number of optical vortices. In other words, if the opti- 
cal Reynolds number indeed corresponds to its fluid counter- 
part, then the quantization of the optical vortices would play 
an analogous role to fluid viscosity. This view seems to be 
echoed by other researchers in the field of superfluids,^^^' 
although we must stress that it is still an open problem as to 
what extent the quantization effect resembles the viscous ef- 
fect.4i 

A. Zero-flux boundary conditions, boundary layers, and 
boundary layer separation 

In classical fluid dynamics the "no-slip" boundary condition 
is most commonly used, and restricts the total velocity and 
hence the total flux to be zero at the boundary. For fluid flow 
above a surface, the velocity shear introduced must be bal- 
anced by a viscous stress, resulting in a boundary layer that 
connects the zero velocity at the boundary to the flow veloc- 
ity above the boundary in an asymptotic expansion."*^ For the 
nonlinear Schrodinger equation, the boundary condition of an 
impenetrable object can be specified by a very low refractive 
index region, which restricts the optical intensity to be zero at 
the surface'*^ due to total internal reflection. Even though the 
tangential velocity can have a non-zero value at the surface, 
both the normal and tangential components of the flux must 
be zero there. This can hence be viewed as a zero-flux "no- 
slip" boundary condition. An optical boundary layer analo- 
gous to the viscous boundary layer in classical fluid dynamics 
is also formed."*^ See Fig.|3lfor a graphical illustration of the 
similarities between a viscous boundary layer and an optical 
boundary layer. 
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Fig. 3. Comparison between a viscous boundary layer and an 
optical boundary layer. 

For a viscous fluid flow past an obstacle, as the Reynolds 
number increases, the boundary layer begins to separate and 
vorticity is convected behind the obstacle. An analogy in the 
dynamics of the nonlinear Schrodinger equation, in the form 
of vortex nucleation on the boundary, is also predicted,^^ 
and in the case of large objects, the instability of the optical 
boundary layer also depends on the optical Reynolds num- 
ber ^ defined in Eq. ( 140^ '^'' much like the viscous boundary 
layer separation. 

B. Dissipation of eddies 

Another important effect of viscosity is the dissipation of 
small-scale structures in turbulence. An analogous effect in 
nonlinear Schrodinger equation is the emission of sound 
waves two vortices are close to each other^"" and the gener- 
ation of Kelvin waves in the process of vortex line recon- 



nections.^^ The radiation of acoustic energy in both cases 
must cause a damping of the high-spatial-frequency convec- 
tion within the optical beam, and the effective Reynolds num- 
ber is again estimated to be equal to the optical Reynolds 
number.^^'^^'^' 



C. Kdrmdn vortex street 

The Karman vortex street is a famous viscous fluid phe- 
nomenon, in which alternate fluid vortices are emitted from 
the back of an obstacle to the flow of a viscous fluid, when the 
Reynolds number increases beyond a certain threshold. ^^'^^ 
Using the numerical vortex blob method, Chorin first simu- 
lated such a phenomenon for a cylinder obstacle and obtained 
good agreement with experimental data.''^*' Since an optical 
beam diffracting past a low refractive index region would also 
emit optical vortices and the vortices interact like vortex blobs 
in a self-defocusing medium, we performed a numerical ex- 
periment of the nonlinear Schrodinger equation to investigate 
if we would observe a similar phenomenon for nonlinear op- 
tics. 
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Fig. 4. Setup of numerical experiment (not-to-scale). 

The numerical setup is sketched in Fig. |3 A big opti- 
cal beam is assumed to propagate at an angle to an ellip- 
soid cylinder, with a refractive index much lower than the 
surroundings to act as an impenetrable object, in a self- 
defocusing medium. The length of the long axis of the el- 
lipsoid cross section is assumed to be W, and the short axis is 
assumed to be one-fifth of W throughout the simulations. The 
two-dimensional nonlinear Schrodinger equation is solved 
using the Fourier split-step method,^'' which implies a peri- 
odic boundary condition for the optical beam. This should not 
affect the qualitative behavior of the dynamics, if the optical 
beam is much bigger than the object. In all of the simula- 
tions, the Mach number ./^ is fixed at 0.4, while the optical 
Reynolds number M is varied. Figure|5lplots the intensity of 
the optical beam at a normalized propagation distance i^ = 10 
for an optical Reynolds number.^ = KW = 12.8. Optical vor- 
tex solitons are created on the top and bottom side of the low- 
refractive-index region, and they interact in such a way that 
resembles the phenomenon of twin vortices behind an obsta- 
cle in a low-Reynolds-number viscous fluid flow. 

Figurel^plots the flux J = (v/*Vi//— \(/'V\(/*)/2i and Fig.El 
plots the momentum vorticity V x J. Both plots confirm the 





Fig . 5 . The intensity of the optical beam at a normalized prop- 
agation distance C = 10, for ^ = 0.4 and ^= 12.8. The dark 
ellipse is the low-refractive-index region that acts as an im- 
penetrable object. Optical vortex solitons are seen to be cre- 
ated on the top and bottom side of the ellipse. While the con- 
vection of the solitons behind the object resembles the twin 
vortices behind an obstacle in a viscous fluid flow. 

similarity between the numerically observed dynamics and 
the phenomenon of twin vortices in a viscous fluid flow. 




Fig. 6. A vector plot of the flux J at ^ == 10, for .-# = 0.4 
and S$, ~ 12.8, which confirms the similarity between the nu- 
merically observed dynamics and the phenomenon of twin 
vortices. 



Figures|8l|9| andllOlplot the intensity, flux, and momentum 
vorticity of the optical beam respectively, at a longer propaga- 
tion distance i^ = 20 for the same parameters. The qualitative 
dynamical behavior of vortices staying behind the object is 
essentially unchanged. 

We now raise the Reynolds number \a M = 25.6 and per- 
form the numerical experiment again. As seen from Figs. II II 
[T2I and[0] the optical vortex solitons become smaller and 
more abundant, but at ^ = 10 the phenomenon of twin vor- 
tices behind an obstacle is again observed. 

At ^ = 20, however, significant instability in the twin vor- 
tices develops, such that the spatial symmetry between the 
upper plane and the lower plane is broken, and alternative 
bunches of optical vortices begin to be emitted from the back 
of the object. Figures [3] El and[^plot the intensity, flux 
and vorticity at ^ = 20 respectively, which demonstrate a be- 
havior strongly resembling the famous Karman vortex street 
phenomenon. 

Due to computing power constraints, we are only able to 
simulate low-Reynolds-number flows, which we do not ex- 
pect to quantitatively reproduce viscous fluid dynamics. We 



Fig. 7. A plot of the momentum vorticity V x J at ^ = 10, for 
Ji = 0.4 and ^ = 12.8. A white dot indicates that flie vor- 
tex has a positive topological charge and a black dot indicates 
that the vortex has a negative charge. The plot shows the sim- 
ilarity between the numerically observed dynamics and the 
phenomenon of twin vortices. 




Fig. 8. The intensity of the optical beam at a normalized prop- 
agation distance C = 20, for Ji = 0.4 and .^ = 12.8. The 
qualitative dynamical behavior is essentially unchanged from 
that shown in Fig.|5] 







Fig. 9. A vector plot of the flux J at ^ = 20, for ^ = 0.4 and 

^=:12.8. 




Fig. 10. A plot of the momentum vorticity V x J at ^ = 20, 
for ^ = 0.4 and ^=12.8. 




Fig. 11. The optical intensity at ^ = 10, for ^ = 0.4 and 
■^ = 25.6. The vortex solitons are observed to be smaller, and 
the phenomenon of twin vortices is again observed. 
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Fig. 12. The flux J at C = 10, for J( = 0.4 and .^ = 25.6. 




Fig. 13. The momentum vorticity V x J at ^ 
0.4 and ^ = 25.6. 
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Fig. 14. Optical intensity at ^ = 20, for Ji = 0.4 and ^ = 
25.6. The twin vortices become unstable and detach alterna- 
tively from the object. 
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Fig. 15. Flux at C = 20, for Ji ^ 0.4 and M = 25.6, which 
shows a flow pattern strongly resembling the Karman vortex 
street. 




Fig. 16. Vorticity at C = 20, for .^ = 0.4 and .^ == 25.6, 
which confirms that the alternate bunches of vortices indeed 
have the right charges that resemble the Karman vortex street 
phenomenon. 



have to use an ellipsoid cylinder in the numerical experi- 
ments, instead of the more conventional circular cylinder, 
to artificially generate more optical vortices, and the Mach 
number is a little too high for compressional waves not to 
play a significant role in the dynamics. With all that said, us- 
ing the nonlinear Schrodinger equation, we are still able to 
qualitatively demonstrate, for the first time to our knowledge, 
two well-known viscous fluid phenomena, namely, the forma- 
tion of twin vortices behind an obstacle, and the symmetry- 
breaking instability of the twin vortices that leads to the 
Karman vortex street when the Reynolds number is increased. 
Compared with previous claims of observing the Karman 
vortex street in nonlinear optics numerically^ or experimen- 
tally,'^'^ our numerical results demonstrate an unprecedented 
level of correspondence between nonlinear optical dynamics 
and the Karman vortex street phenomenon, thanks to the pres- 
ence of a much larger number of optical vortices in our simu- 
lations. 

D. Kolmogorov turbulence 

The striking similarities between nonlinear optics and vis- 
cous fluid dynamics are not limited to low-Reynolds-number 
two-dimensional problems. As the Reynolds number is fur- 
ther increased to the order of a million, the viscous fluid 
flow enters a turbulent regime. Since this regime is highly 
chaotic, only statistical signatures can be reproduced in a tur- 
bulent fluid flow. A well-known signature of turbulence is the 



Kolmogorov energy spectrum,^^ derived under the assump- 
tion that a "steady state" is reached when the macroscopic- 
scale fluid flow continuously generate finer spatial structures 
via convection and viscosity dissipates the smallest struc- 
tures. As viscosity plays a significant role in the Kolmogorov 
turbulence spectrum, it is surprising to see that numerical 
simulations of the three-dimensional nonlinear Schrodinger 
equation also reproduce the Kolmogorov spectrum at high 
Reynolds numbers, and the vorticity dynamics of the "su- 
perflow" described by the nonlinear Schrodinger equation re- 
sembles that of the viscous flow, in which vortex reconnection 
events play a major role. ^^'^^ 

The dissipation of the smallest spatial structures in a su- 
perflow is speculated to be the Kelvin waves produced by 
the natural motion and reconnections of vortex filaments,^''^''^ 
and the corresponding Reynolds number is again speculated 
to be Si ~ KW f''^ Numerical and theoretical analyses of the 
so-called "quantum turbulence" exhibited by the nonlinear 
Schrodinger equation all reveal striking similarities between 
quantum and classical fluids, and it is argued that the study of 
quantum turbulence could lead to a better understanding of 
turbulence in normal fluids.'**^ 

4. The split-step method 

While a nonlinear optical system shows promise for com- 
puting Euler and Navier-Stokes fluid dynamics, it also poses 
serious technical challenges. Ideally one would like to have 
a configurable nonlinear material with low loss, anomalous 
group-velocity dispersion, high defocusing nonlinearity, and 
three-dimensional co-propagating boundaries. One may only 
be able to find separate materials or optical devices, each of 
which performs only some of the functions. Moreover, par- 
asitic effects such as loss, two-photon absorption, and high- 
order dispersion can be detrimental to the accuracy. To com- 
bine different devices and periodically compensate for para- 
sitic effects, we hereby propose the "split-step" method, the 
inspiration of which comes from the numerical Fourier "split- 
step" method.^" Consider the general nonlinear Schrodinger 
equation. 



the error order, 



d\\l 



N 



(41) 



n=l 



where propagation effects and boundary conditions are ex- 
pressed in terms of operators H„. The formal solution is 

V/(C+AC)=exp(|^^^^f H„dC')r(C)- (42) 

But if A ^ is much smaller than 1 /H where H is the magnitude 
of the operators, by virtue of the Baker-Hausdorff formula we 
have 

V/(C + AC) = n^xp (^»AC) r(C) + 0{h'^aC^). (43) 

Each of the propagation effects can hence be applied sepa- 
rately to an optical pulse, with a quadratic error term. A sym- 
metrized version of the spht-step method can further reduce 



V/(C + AC) = n exp(//,„^) nexp(//„^)r(C) 

m=N ft= 1 



o{wm;') 



(44) 



The split-step method is not unlike the proof of a quan- 
tum computer being able to simulate any quantum systems. ^^ 
Whereas it is difficult to find a quantum device that performs 
the exact Hamiltonian of the quantum system of interest, it is 
possible to approximate the Hamiltonian in small time slices. 
Similarly, in a metaphoric optical computer, one can form a 
unit cell of a "meta-material" by combining a slice of defo- 
cusing material, a slice of material with anomalous group- 
velocity dispersion, a slice of ultrafast phase modulator to 
apply the three-dimensional boundary conditions, and a gain 
medium to compensate for loss. The optical beam can loop 
through the unit cell multiple times in a cavity, so that the 
outcome will approximate the true solution as if we had an 
ideal medium. See Fig.[n]for a graphical illustration of the 
method. 
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Fig. 17. Sketch of a split-step optical system that approxi- 
mates the ideal nonlinear medium. 

The split-step method has the additional advantages that 
each subsystem can be tunable and easily substituted with an- 
other material or device, and the pulse evolution can be mon- 
itored more easily. The magnitude of each effect can be tuned 
by simply changing the propagation length in each device. 
In exchange of configurability we have sacrificed some ac- 
curacy due to discretization errors and instability. The com- 
putation speed may also be reduced by a large but constant 
fraction, as the pulse may spend most of its time on simply 
propagating from one device to the next and not performing 
the core computation by nonlinear propagation. The spUt-step 
method, however, does not detract from the inherent paral- 
lelism in the computation, as the transverse dimensions are 
not discretized. 

5. Conclusion 

In conclusion, we have used a variety of theoretical and nu- 
merical methods to show that self-defocusing optical propa- 



10 



gation has a converging correspondence with Euler fluid dy- 
namics and a striking similarity with Navier- Stokes fluid dy- 
namics. We have numericafly shown that the interactions of 
a large number of optical vortex solitons are able to simu- 
late two well-known viscous fluid phenomena. We have also 
proposed the split-step method, a way of practically imple- 
menting the metaphoric optical computer. 

There are serious technical challenges if a metaphoric com- 
puter is to become useful for computing fluid dynamics, es- 
pecially three-dimensional fluid dynamics problems, as tech- 
niques for the complete specification and characterization of 
the spatiotemporal optical field are still in their infancy. The 
speed, configurability, and parallelism of a metaphoric optical 
simulator nonetheless promise vast advantages over conven- 
tional numerical simulations. 

Since photons are quantum objects, optical propagation 
would also inherently compute the quantum dynamics of 
bosons, and may therefore be used as a metaphoric simula- 
tor of quantum fluids, such as superfluids, superconductors, 
and Bose-Einstein condensates. In this way the advantages of 
a metaphoric computer and those of a quantum computer are 
combined, and only then the classical and quantum comput- 
ing capabilities offered by photons would truly be exhausted. 
This extension of metaphoric optical computing will be a sub- 
ject of future work. 
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